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In view of the recent WA98 data of tv° spectra from central 
Pb + Pb collisions at the CERN SPS, we analyze the produc¬ 
tion of neutral pions for A+A collisions across the periodic ta¬ 
ble at y/s = 17 A GeV and 200 A GeV within the framework 
of the parton-cascade model for relativistic heavy ion colli¬ 
sions. The multiplicity of the pions (having p± > 0.5 GeV/c) 
in the central rapidity region, is seen to scale as ~ ( N par t )“, 
where N par t is the number of participating nucleons, which 
we have approximated as 2 A for central collisions of identical 
nuclei. We argue that the deviation of a (~ 1.2) from unity 
may have its origin in the multiple scattering suffered by the 
partons. We also find that the constant of proportionality in 
the above scaling relation increases substantially in going from 
SPS to RHIC energies. This would imply that the (semi)hard 
partonic activity becomes a much cleaner signal above the 
soft particle production at the higher energy of RHIC, and 
thus much less dependent on the (lack of) understanding of 
the underlying soft physics background. 


Recently, the WA98 Collaboration has published fj]] 
data for the production of neutral pions up to transverse 
momenta of p± ~ 4 GeV/c, in central Pb + Pb collisions 
at 160 A GeV/c incident momentum, corresponding to 
y/s ~ 17 A GeV. Two most interesting features of these 
data emerge when compared to corresponding data from 
pp collisions and collisions involving lighter nuclei §: a) 
an approximate invariance of the spectral shapes, i.e., a 
near indepence of the slope of the neutral pion p± spec¬ 
tra; b) a simple scaling of the n° with the number of 
participating nucleons, if the number of participants is 
large ( > 30). 

Soon after its publication the (preliminary) WA98 
data |'J], was analyzed in two interesting papers, one 
by Wang || and the other by Gyulassy and Levai [||, 
with regard to the contribution to n° production at 
larger pj_>,2 GeV/c due to (semi)hard parton scat¬ 
terings and associated minijet production. Astonish¬ 
ingly, Wang was able to reproduce the WA98 data for 
the invariant 7r° cross-section by simply using the per¬ 
turbative QCD cross-section for a single parton-parton 
scattering convoluted with the nuclear overlap func¬ 
tion and a quark/hadron fragmentation function, and a 
parametrized accounting of the pp -kick suffered by the 
partons. On the other hand, Gyulassy and Levai made a 
detailed numerical simulation using the event generator 
HIJING, and claimed that the consideration of multiple 


scattering effects of partons are crucial and not at all in 
conflict with the WA98 data. 

In the present paper, we analyze the WA98 data from 
yet another angle [| | . We would like to demonstrate that 
the observed scaling of the pion yield in A + A collisions 
as well as the insensitivity of the p± slope to a variation 
of A, does not at all imply the absence of dense matter 
effects on the particles, nor does it necessarily lead to 
the conclusion of Ref. 0 that the produced matter is 
thermally equilibrated. Gyulassy and Levai have already 
hinted this, but our findings go even further. In fact, a 
significant partonic cascade activity which scales roughly 
with A 4 / 3 (rather than with A) does not lead to a con¬ 
tradiction with the WA98 data, rather, we believe that 
it is confirmed by the WA98 results. 

In order to elucidate this apparent mystery, we use the 
event generator VNI Q which embodies the physics of 
the parton-cascade model 0 for ultra-relativistic heavy- 
ion collisions. The model attempts to describe the nu¬ 
clear dynamics on the microscopic level of particle trans¬ 
port and interactions, by evolving the multi-particle sys¬ 
tem in space-time from the instant of nuclear overlap all 
the way to the final-state hadron yield. For details we 
refer the interested reader to Refs. [||. Here it suffices to 
summarize the essential elements embodied in the model, 
namely: (i) the initial minijet production through liber¬ 
ation of partons from the colliding nuclei by means of 
(semi)hard parton scatterings; (ii) the subsequent parton 
cascading consisting of multiple gluon emission and suc¬ 
cessive rescatterings; (iii) the coalescence of final-state 
partons to color-neutral pre-hadronic clusters; (iv) the 
decay of clusters into primary hadrons; (v) the subse¬ 
quent hadron cascading consisting of reinteractions of 
hadrons and clusters with resonance production and de¬ 
cay. 

In this model, ‘dilute’ collision systems (involving 
beams of protons or light nuclei) naturally show a very 
different dynamical evolution than ‘dense’ collision sys¬ 
tems (head-on collisions of heavy nuclei). Whereas in 
the former case multiple particle interactions are absent 
or negligible, in the case of heavy ions, the space-time de¬ 
velopment is characterized by multiple scatterings of par- 
tons, and later also of hadrons which are formed when 
the partonic system hadronizes. Moreover, the parton 
cascading during the early stage has rather distinct fea¬ 
tures, as compared to the hadron cascading during the 
late stage of the collisions. Multiple interactions among 
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produced hadrons mostly alter only the momentum dis¬ 
tributions and do not lead to a substantial increase 
of particle production. The interactions among partons, 
on the other hand, have a very different effect on par¬ 
ticle production. Perturbative QCD tells us that pri¬ 
mary (semi)hard scatterings of partons with momentum 
transfer q± > 1 — 2 GeV/c contribute the bulk of minijets. 
Moreover, the so kicked-out partons are excited off-shell 
by an amount Q 2 ~ which they tend to shake off 
by gluon brems-strahlung. Thus, even in the absence of 
partonic rescatterings, the additional emission of gluons 
leads to an increase of the parton multiplicity on top of 
the number of those quarks and gluons which are liber¬ 
ated by the primary scatterings. In the parton-cascade 
model the characteristics of parton multiplication is fur¬ 
ther pronounced due to secondary interactions: Firstly, 
the initially kicked-out primary partons may rescatter 
and receive additional momentum transfer that again 
feeds the gluon emission. Secondly, all the newly pro¬ 
duced off-spring of secondary partons increase the local 
density and can themselves re-interact. However, upon 
hadronization, the partonic color charges have to reor¬ 
ganize themselves to color-singlet composite objects that 
are the seeds of the emerging hadrons. Clearly, in or¬ 
der to make a hadron of mass rrih , the total invariant 
mass of the recombined partons must be at least equal 
to uiji , and consequently it requires in the mean more 
than just two partons to satisfy this kinematic condition, 
especially since the bulk of gluons from radiative emis¬ 
sion piles up at the low-energy end of the perturbative 
regime. Nevertheless one would expect a proportionality 
between produced partons and resulting hadrons. 

Indeed, the cluster-hadronization scheme |J, employed 
here to convert the outcome of the parton cascade into 
hadronic states, provides that the number of produced 
pre-hadronic clusters, and hence hadrons is proportional 
to the number of partons present in the system. Thus, 
the production of particles should reflect the extent of 
partonic scattering throughout the nuclear collision. This 
can be a very important consideration as the number of 
collisions among the primary partons has been shown to 
scale linearly with the number of the participating nucle¬ 
ons in the nuclear collision [[lo| . Thus a deviation of the 
number of produced hadrons from a similar linear scal¬ 
ing would be a direct confirmation of multiple scatterings 
taking place in the wake of relativistic collision of nuclei. 

A simple consideration may illustrate these features of 
particle production. Let x denote the number of partons 
in each nucleon, and let each parton suffer u collisions 
during the partonic stage. Assuming that each virtual 
parton radiates r partons, we see that the number of 
produced partons will vary as 

Apartons OC V (1 + r) X A , (1) 

and if, as proclaimed, 


Ahadrons x Ap ar tons , (2) 

one realizes immediately that if the partons interact only 
once, the multiplicity of the partons, and hence the mul¬ 
tiplicity of hadrons, will scale as A. It is also clear that if 
every parton interacts with with every other parton then 
v oc A, and the number of materialized partons would 
scale as A 2 . That can happen, if the system would live 
for an infinitely long time. However, this is not the case. 
Rather than that, in relativistic heavy ion collisions, the 
partonic matter will expand, dilute, and eventually con¬ 
vert into hadrons. Thus a given parton may undergo 
v ~ R/X interactions; where R is the transverse size 
of the system and A is the mean free path of the par- 
ton. Noting that R ~ A 1 / 3 , we immediately see that the 
number of materialized partons, and hence the number 
of produced particles would scale as ~ A 4 / 3 . An exper¬ 
imental verification of this scaling behaviour could be a 
direct manifestation the formation of a dense partonic 
matter! 

We shall demonstrate now that these simple consid¬ 
erations are indeed confirmed by a detailed simulation 
with the event generator VNI on the basis of the parton- 
cascade/cluster-hadronization model. We first consider 
the recently measured transverse momentum distribution 
of 7r°-production in central collsions of Pb + Pb at CERN 
SPS obtained by the WA98 collaboration jj|. To make 
contact with the experimental data, the simulations were 
done for the range of impact parameters 0 < b < 4.5 
fm, which corresponds to 10% of minimum-bias cross- 
section. The result of our model calculation, shown as 
the solid histogram in Fig. la, is seen to be in decent 
agreement with the experimental measurements. The 
model results do not include the final-state interaction 
among produced hadrons yet, but it is likely || that the 
agreement will further improve once the effect of cascad¬ 
ing hadrons is included. The dashed histogram in Fig. la 
corresponds to pp collisions at y/s = 17 GeV, and scaled 
accordingly. Comparing Pb+Pb with pp, clearly exhibits 
the well-known enhancement of the production of pions 
having large p± due to multiple scatterings among the 
partons ||10| . 

As mentioned already, one of the interesting observa¬ 
tions of the WA98 results is near independence of the 
slope of the p± spectra for the pions once the number 
of participants is large. In Fig. lb we plot our results 
for the p± spectra of 7r°’s for a number of central A A 
collisions at y/s = 17 A GeV for various A + A systems 
from A = 16 to A = 197. One observes that they are al¬ 
most identical in shape with a universal slope for p± ,<1.5 
GeV/c. On the other hand, the deviations appearing at 
larger p± for heavier systems are indicative of enhanced 
multiple scattering there. Similar results (not displayed 
here) were obtained at RHIC energies. 
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7V,o cx (N part ) a , (3) 

where N part = 2 A is the number of participating nucle¬ 
ons, and a being extracted as: 

f 1.16 at i/s = 17 A GeV ... 

a ~ \ 1.23 at yfs = 200 A GeV ' 14j 

In view of the uncertainties, we may say that a ~ 1.2 is 
a fair number that characterizes our model results. This 
is quite similar jn| to the observation of the WA98 ex¬ 
periment at y/s ~ 17 A GeV, namely [jf] a ~ 1.3, with 
the number of participating nucleons estimated from the 
impact parameter, in collisions involving lead nuclei. 
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FIG. 1. a) Top: Transverse mass spectra of neutral pi- 
ons in central (10% of minimum bias) collision of 158 AGeV 
Pb + Pb collisions. The solid histogram represents our result 
from the parton-cascade/cluster-hadronization model. The 
dashed histogram corresponds to pp collsions at E crn = 17 
GeV, scaled accordingly, b) Bottom: Transverse momentum 
spectra of neutral pions in central collisions of identical nuclei 
at E crn = 17 A GeV. The symbols are our results for various 
systems, from O+O up to Au+Au. All results are normalized 
to the case A = 16. 

In order to verify this scaling more closely, we have 
calculated, as a function of the nuclear mass number 
A, the production of 7r°’s in the central rapidity region 
(—0.5 < y < 0.5) having transverse momenta p± > 0.5 
GeV/c. The latter choice minimizes the influence of pi¬ 
ons having their origin in decay of resonances. This kine¬ 
matic window was motivated p ] by the WA98 collabora¬ 
tion in their measurement of the n° yield. Fig. 2a dis¬ 
plays the simulation results for central A + A collisions 
across the periodic table, at CERN SPS center-of-mass 
energy yfs = 17 A GeV, while Fig. 2b shows the same 
for RHIC energy y/s = 200 A GeV. The solid lines are 
fits to the model results, represented by the symbols, and 
scale as 
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FIG. 2. a) Top: Mass number scaling of pion production 
in the central rapidity region at CERN SPS energies. The 
symbols represent the results of our simulations, and the solid 
curve is a fit to these results, b) Bottom: Mass number scaling 
of particle production in the central rapidity region at BNL 
RHIC energies. The symbols represent the results of our sim¬ 
ulations, and the solid curve is a fit to these results. Dashed 
lines would correspond to a linear scaling, dotted-dashed lines 
to a scaling with the nuclear overlap Xaa(5 = 0). 
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For comparison, the dashed lines correspond to a lin¬ 
ear scaling ~ ( N part j 1 ' 0 , whereas the dashed-dotted lines 
indicate a hypothetical scaling with the nuclear over¬ 
lap factor Taa{ b = 0) ~ (Apart) 1 ' 42 - A linear scal¬ 
ing would reflect a single-collision situation, and a scal¬ 
ing with Taa would indicate a Glauber-type multiple- 
collision scenario, in which nucleons suffer several colli¬ 
sions along their incident straight-line trajectory, with¬ 
out deflection but with energy loss. Comparing the three 
curves, we can conclude that our simulation results rise 
significantly slower than with Taai because firstly, the 
particles change direction through the collisions, and sec¬ 
ondly, they are subject to a collision time of the order of 
the inverse momentum transfer, during which they can¬ 
not rescatter. On the other hand, our calculated n 0 yields 
grow much faster than linear with A, due to multiple scat¬ 
terings. Hence, the scaling of the 7 r° yield in the model 
may be interpreted as a collision-meter that indicates the 
extent of multiple interactions on the parton level (recall, 
that we did not include hadronic final-state interactions). 
Such an interpretation would also explain the (slightly) 
larger value of a, and also the constant of proportionality 
(cf. Fig. 2a and 2b) i.e. the stronger increase, at RHIC 
energy as compared to the results for CERN SPS energy. 
In the former case one would expect that multiple col¬ 
lisions are further enhanced due to a higher density of 
produced particles. This, we find, results in a increase 
larger than a factor of 2 in the number of 7r°’s between 
CERN SPS and RHIC. 

In summary, we have demonstrated here that the ob¬ 
served scaling with A a with a — 1.3 of the number of 
produced particles in heavy-ion A + A collisions, as well 
as the approximate shape-independence of the transverse 
momentum spectra, are satisfactorily reproduced by the 
parton-cascade / cluster-hadronization model, though 
with a slightly smaller value of a ~ 1.2. The quark- 
gluon multiplicity from the early stage of the collisions 
due to initial minijet production plus subsequent cascad¬ 
ing with multiple scatterings and gluon emissions, leads 
to a final-state particle multiplicity which approximately 
scales with the number of nucleon collisions (i.e., ~ A 4 / 3 ), 
as inferred by the WA98 experiment using high p± 7r°’s. 

The key result of our model simulations is that multiple 
parton scatterings contribute significantly in collisions in¬ 
volving heavy nuclei. We have checked that the initial 
minijet production component alone, without subsequent 
cascading and reinteractions of the minijets, would yield 
for heavy ions a 7r° yield as for the pp case in Fig. la, 
which would clearly underestimate the WA98 Pb + Pb 
data. 
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